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ABSTRACT 


The reversion-of-series method is extended to the 
s - domain by utilizing non-linear Laplace transforms. 
The reversion of series in the s —- domain is applied to a 
non-linear differential equation and approximate solutions 
are obtained. The approximate solution is modified for 
the case where the steady state is a constant value by 
calculating the exact steady-state value and applying it 
to the reversion approximation. The non-linear differen- 
tial equation considered is Duffing's equation with a 
danpingwcerm and sinusoidal and constant forcing functicnee 


The theoretical solutions are compared to machine solutions. 
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IT. INPRODUCTION 


Utilization of the Laplace transformation to obtain an 
accurate analytical solution of a non-linear differential 
equation has been prevented, in past, by the inability to 
obtain an appropriate expression for the Laplace transform 
of the non-linear terms. Baycura [Ref. 1] developed an 
expression for the Laplace transform of certain non-linear 
terms and Brady [Ref. 2] obtained a formulation for the 
general case. The non-linear transform is used in con- 
junction with the reversion-of-series method to obtain an 
approximate solution for a non-linear differential 
equation. 

The approximate solution is formulated and examined 
in detail by comparison to a machine solution. A solution 
obtained by modifying the approximate solution is also 


derived and examined in detail. 


TI. TRANSFORM EXPRESSION FOR INTEGRAL POWERS OF A FUNCTION 


The general formula for the Laplace tse tormeeiaa 


function raised to an integral power [Ref. 3] is: 


fic es. ere) | GL) 





Where Nis Aspesi tive minterer treater Ieee Zener 


TTil. REVERSION OF Seen 
If a series is represented by 
= ay + a i +a y> + a my + a y? + (a. # 0) (2) 
1 2 3 4 5 eet 1 


the coefficients of the series 


P 2 3 A 5 
y = A,x + A,X + A 2x + Ay x + Agx Ecos (3) 
are given by [Ref. 4]: 
if q 1 
A, =— A, =--—-> A, = —= (2a - aa.) 
al a5 2 a8 3 a7 o dls 


_ 2 2 3 
Ay = 7 (58858. - a) a, - Dan ) 
a 
1 
— @ o ee 4 5 2 
Ag = ~ 9 (6a, Ana) + 3a, a. + lia, - 4,785 - 2la,a, 2.2) 
di 


Utilizing the relationship expressed in Equation (1) the 


Laplace transform of Equation (2) becomes: 


eS) = a,Y¥(s) tassY°(s) +a,5°¥3(s) +a,s°¥ ‘(s) ee Gi) 
Let: 

_ » _ 2 - 3 
Dy = a4 Ds, = 58 D3 = 228 Dy = 2,8 


Equation (4) can then be rewritten as: 


X(s) = b,Y(s) + boY°(s) + baY9(s) + byY(s) oy G5 ) 


The Laplace transform of Equation (3) is: 


Y(s) = A.X(s) + A sX“(s) + A sere > (ia) + A Bee le) 5 (6) 
1 2 3 4* 
Let: 
B, = @2 Be = Ans Be = fh =" Bae = A 52 
iL 1 2 e 3 = 4 mM 


Equation (6) can then be rewritten as: 
vig) 8 E.G) > Bocls) Hyon Ene le (7) 
ci 2 2 \ om « 


The coefficients of Equation (7) can be expressed in terms 


of the coefficients of Equation (5): 


dl 2 di 2 
CS ba. = = => Bi, = — ee (2b Fie Win) 
eee Sl Bee 2 ee LS 
dL i 
at 2 3 
By = it (5b, bb 2 - bo Dy = DD. ) 
1 
Pt dll 2 eae 4 5 Z 
Be = coe bob) i 3b, D3 1 14b,, = DS Wis = elb,b,, bd.) 
i 


Thus, Y(s) can be expressed as a function of X(s) by the 


reversion of Equation (5) into Equation (7). 


IV. GENERAL SOLUTION OF A NON-LINEAR DIFFERENTIAL EQUATION 


The technique illustrated in the previous section will 


now be used to solve Duffing's equation with a damping 


term and initial conditions. The equation is expressed as: 
ay dy : 3 

Ae eae * Og yee =e) (8) 
t 


where a, we and h are constants and x(t) is a forcing func-— 


tion. The Laplace transform of Equation (8) is: 


"EG + as + we )¥(s) ~ ns“¥7(s) = X(s) +sy(0) + Z| +ay(0) (9) 
0 


Let: 

Pees 2 
by =s + as + Wo 
D3 = hs¢ 


Z(s) = X(s) + y(0)(sto) + 2% 
dat 0 
Equation (9) can then be written as: 


b,Y(s) + baY°(s) = IGS. (10) 


Consider Equation (10) as series with oe equal to zero for 
all n except n equal to one and n equal to three. 


By reversion of series: 


vis) = —-2(s)_  _ _ns®z3(s)_, 3n®s'z%(s)_, — 
a eee 
(s“+os+ws) (s“tastwe) (s“+astue) ! 


The denominator of each term in Equation (11) can be 


Pactore Gia hus: 


Z gag oO MR ee eee 
s“ + as t+ Wy = | 8 + 5 + 5 s + 5 ; (125) 
Let: 
i 
cee _ 
= x + 5 (a Aw) Cie 
oe. 
= "a _ 
Equation (11) can then be rewritten as: 
OS AN Is 
Hae Zs) = hs Z~(s) penile) 7 (s) 2 we (15) 


of ae eee 
(umidstB) (sf8) (SHB) an) | guy! 


Equation (15) represents the general solution in the 

s -— domain. Inspection of Equation (15) reveals that the 
form of the time solution depends on the nature of the 
forcing function and the relationship between 0° and hue. 
Further, if the forcing function is Laplace transformable 


and if h is equal to zero, the solution is exact. 


ane 


V. APPLICATION OF THE GENERAL SORUIION 


The general solution just derived will now be applied 
fer zero initial conditions with constant and ‘Saimieeida 
forcing functions. The time solutions are given in terms 
of numbered coefficients. The formula for each coefficient 


is listed in the Appendix. 


A. CONSTANT FORCING FUNCTION 

The inverse Laplace transform of Equation (15) yields 
three different solutions depending on the nature of the 
Radical in Equations (13) and (14). 


1. Radical Equal to Zero 


For the radical to be equal to zero, then: 


2 
am = dy (16) 
0 
and 
= ES 
= = Ss Cie) 
The general solution in the s - domain is then: 
3 29) 
Ine Site Cc 
Y(s) = ——— - +t te. (18) 
s(s+A)*° Ss Gst) s(sta)* 


Where Carsmune constant Forcing at uneraen. 
The inverse Laplace transform of the first three terms of 


Equation (18) yields: 


de. 


7 Se 
y(t) =K, + (K5+K.t Je + Ky + (Ke+Ket+k 


2 13s 2 he 
7b arta Ki gt Je 


(129 


where the Kas i = 1, 2, 3, .+., are defined in the Append 


An examination of Equation (19) indicates that the initial 
conditions are satisfied. The effect of the various param- 
eters on the solution can be seen by examining the steady- 


state solution which can be expressed as: 











2 2 7 ae 
<a sale he “ he 
a sheers eye] ewes p pS 
0 0 0 0 


It should be noted that only the first three terms of 
Equation (20) are utilized in Equation (19). Since Equation 
(20) is an alternating infinite series, a condition for 
convergence [Ref. 5] is that the magnitude of each successive 
term must decrease. This implies that the applicability 

of the reversion method is limited in the steady state to 


the situation where hee is small compared to (we)?. 


Using the arbitrary values c = 1.0, ao = 4.0, we = hea 
and h = 1.0, the solution represented by Equation (19) 
is compared in Table I to a machine solution which was 
obtained using Runge-Kutta Adams - Moulton with error check. 
As the results indicate, a relatively accurate approximate 
SOuuctOn was Obtained. in Table Ii the solutions are cone. 
pared for the coefficient of the non-linear term equal to 


ten. In this case the solution was a good approximation, 


but not as accurate as in the previous example. Examination 


ale 


TABLE I 


w=4.0, we=l 0, h=1.0., cadno 





asi Linear Réeversaon Machine EpRer 
Choe O°. OnGrele Coe 0.0000 .0000 
Che 0.0661 0661 Shs liou.o) .0000 
EO 0.1485 ~1485 0.1484 C01 
digs Oreuve Teele 0.1997 .0005 
2m Ome oc (al 2269 0.2259  COH0 
Zep 0.2399 2394 Ol 2) SS Sige 
50 Oe y5 7 2447 0.2430 Oa: 7 
Sas) 0.2482 2466 Ore oer Ow > 
4.0 0.2492 eal 0.2458 - O@H 3 
4.5 0.2497 2471 0.2461 Ss Oui IL0 
DO 0.2499 2467 0.2462 . 0005 
5+ 0.2499 ~2466 0.2463 a 9,0).0) 5 
6.0 Ore 700 ~2465 0.2463 .0002 
er OTe S010 ~2464 0.2463 mo) 0001 
i. 0 Ores 00 2no 3 0.2463 .0000 
is) Oe25 010 2463 Oh 2 iars . CIO 0 
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TABLE ii 


a=4.0, wo=4.0, h=10.0, c=1.0 





—- Fr WH W 
UI 


CO) eee WON U1. 
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Ise 


Linear 


A OUOLENO, 
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~1485 
we O02 
wae 
SBI, 
se 457 
2482 
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ne 0 0 
qe 50.0 
5 SONG 
42500 
ao Ole 
we) OO 


m9) 010, 


Reversion 
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1997 
BRCM 
2347 
ne Bal 
Re 34, 
S225 
2246 
mee lg 
yee Os 
mece 
22206 
ee) 
woes 
2246 
2268 
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Or, 


Machine 
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0660 
-LATT 
.1956 
2 eM 
Ri foatoual 
Mee Oe 
ween 
pee | 
Welle 208 
TES) 
Deefad 22 3 
DEES. 
oe 
Pees D 
meee 
ee 
5 ES 


S222 


Error 


.0000 
my 6)1 6/10) 1. 
0008 
OCs 
OCs 
. OUZS 
OMe 
0098 
005s 
. Gem 
0006 
ra) 8).2"(0) 
00a 
» OOEE 
. OCs 
.0008 
. OCs 
0043 
yO Is) 1 


of Equation (20) reveals that for h equal to ten, successive 
terms decrease in magnitude, but not as rapidly as in the 
case h was equal to one. If four terms of Equation (20) 
had been utilized in the solution, a steady-state value of 
0.2178 would have been obtained resulting in a better 
eeproximation. Table JII represents a comparison om tie 
s@futions for the forcing funetion equal to ten. Imes 
Gase the solution limits were exceeded and a large posmuame 
number was obtained for the steady state. If an additional 
term of Equation (20) had been included in the solution, a 
larger negative value would have been obtained for the 
steady state. Good results were realized however, for the 
imjittiad porhiken of the tramsient solution, ingicatime the 
applicability of the reversion method in this region. The 
linear solution represented by the first two terms of 
Equation (19) was also included in the tables to illustrate 
a less accurate but simplified approximate solution. 

2. Rad#eal Greater than Zero 


The relationship between 0° and hwé is expressed as: 


2 @ 


Ol >Awo (21) 
Equation (15) is then expressed as: 

3 25 
Y(s) = ———___ _ - —__"°_+ shime + s—sl@zen 


s(stA)(stB) s(stA) '(stB) a(t eee 


The inverse Laplace transform of the first two terms of 


Equation (22) yields: 


IE) 


TABLE III 


gs 10, weal. ai 0, 611,80 





Time 


\O 


a0, 


Co on tf 


= aaa oe 


co) om Ui 


ar Ol’ © 


> Bee 3| 


Linear 


.0000 
Sei) 
.4780 
. 8434 
mlor7 
4850 
.0020 
P70 
moo 90 
.4570 
4818 
4925 
.4969 
~4988 
4998 
.5000 
.5000 
5 UN0L8 
.5000 


500 


dmaly. 
14. 
15. 
oy 


Reversion 


.0000 
Oe 
.4780 
. 8433 
mangiony 
.4807 
. 9560 
.0727 
.8907 
roa 
. 1488 
.9814 
52,0) ,64l 
.9344 
. 7740 
2 Oe 


9169 
3046 
6965 
9043 
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OF 


ys 


0 


on 


Machine 


0000 


oe 


4770 


ool 


mC te 
4085 
ale) ios, 
Oo 
» 5006 
oo 
3551 | 
TOMAS 
Digs 
5567 
oo 
BO Ee 
5 ee 
Sool 
.5568 
5 leye 


Error 


O% 


V2" 
14. 
le ee 


0000 


-0000 
mG) 6) 7 10, 
NE 2 
. 02e8 
. O fae 
. sue 
4847 
» 3901 
0285 
4079 
Sean 
» 3488 
ee i 
cli 
99ey 
ING) 


3601 
7478 
ao | 
Sue 


y(t) =Ky, + Ke’ + Kone) +K 


a0. it 
eye. Je 


2 
23 + (Ky y + K55¢ + Ket tKy ot 


2 oe BiG 
+ (Kog + Kogt + Kot +K.,t Je : (23) 


The inverse Laplace transform of additional terms of 
Equation (22) results in a steady-state solution given by 
Equation (20). The region in which the solution can be 
applied is then the same as discussed in the previous 
section. Using the values ¢ = 1.0, a = 6.0, w. = 4. 0mema 
h = 1.0, the solution represented by Equation (23) and the 
solution obtained considering the first three terms of 
Equation (23) are compared to a machine solution in Table 
IV. The small errors listed indicate that a highly accurate 
approximation was obtained. Table V illustrates the effect 
of increasing the coefficient of the non-linear term to ten. 
A less accurate approximation was obtained. The effect of 
increasing the magnitude of the foreing function was not 
illustrated; however for the parameter values used in 
Table IV if ec had been increased to ten, a steady-state 
value of -1.4063 would have been obtained indicating that 
the solution limits would have been exceeded. The linear 
solution is seen to be a less accurate approximation. 

3. Radical Imaginary 


The relationship between a* and hws see 


lw? > af . (2H) 





TABLE IV 
u=6.0, wh=4.0, h=1.0, c=1.0 
ame pene Y Re Vie s won Maclin Era 
O70 . 0000 0.0000 .0000 Ova 0:00 
DY 5 0).5) S18) 00533 0535 0.0000 
Las, ois 9 Omar 59 jl dees Ou O'CIO ul 
5 21570 OmM569 SSeS 0...0Gi0a 
20 sigs 0.1864 .1860 0.0004 
25 w2066 0.2064 .2058 0.0006 
3) 8, 2204 OR 2 016 209i 0 20019 
B.5 2298 2292 . i 1. Om O11 
4.0 #2802 Ome 355 .2341 0. n0@ak2 
5) 2436 0.2420 2408 0 ..0:0mez 
6.0 2470 0.2448 2438 0 . O1OMKG 
gue 2486 0.2459 2452 O . OCias, 
8.0 2494 0.2462 2459 0.00 
9,0 2497 0.2463 2460 D OOe 
are!,.0 2499 0.2463 2462 Ovi@iOlOal 
am0 2499 0.2462 2462 0.0000 
2.0 200 0.2462 246 3 -0.0001 
3.30 2500 0.2461 2463 -0.0002 
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TABLE V 





a=6.0, wo=4.0, n=10.0, o=1.0 
ene Lames Reversion Machine Eeror 
oO ~0000 O20010 0 .0000 0.0000 
a5 205 33 0205355 a0 5.5 0. Gee 
0 al Loe OL lanes bal S36 Om00:0 3 
IL EO 0.1568 mls) oy: O #0) 
ae. Pilis ions Oeakcs 7 1821 0.0036 
a, 2066 0.2047 1986 0.0061 
5.0 .2204 C185) 2086 0.0079 
Be) 2298 0.2234 2145 0.0089 
4.0 . 2362 0.2269 R279 0.0090 
4.5 ~2406 Olen. 2196 OMMOOIS 
DU 2436 Ose27s Ry aes AG 0.0068 
Shp. »2465 0.2266 BAUS: 0.0050 
6.0 Pee a) Omee 50 ree OPre C0 
O.5 ~2480 Oe223l eee 0.0009 
f M0 .2486 Oma 5 wee 5 =O FOC 
oro .2494 0. 2mr61. .2224 =O OWS, 3 
9.0 ~2497 Ora wee 5 =. 0068 
10 .o 2499 0... caso 3) =O", O16 
Lee!) 5 ORG C220 2225 -0 . Omg 
aS 0 eo Pe ae) Nz A) —O. CaS 


i 


5 (o%—4we) = 

Then 

A= 5+ 5D C25 ) 
and 

B= 5- gD. (26) 


Equation (15) is then expressed as: 


te, eG 
s(stA)(st+B) s(stA) (stB) s(stA)'(s+B) / 


came 
The inverse Laplace transform of the first two terms of 


Equation (27) yields: 


y(t) =Kyg t Kye + Kyoe + Kyg t (Kyy + Kyst + Kypt+Kyotrye 


fo RL +o ene (28) 


+ (Kip +Kygt +Keg 51 


Since A and B are complex, the coefficients of the time- 
dependent portion of Equation (28) are complex. Equation 
(28) can then be rewritten as: 

_ : -At : -Bt 
act ) Kg + (M,tIM,)e + (M,-JjM,)e + Kia 


: P : C , So ae 
+ [(M24+5My) + (Mo+jMg)t + (Mo+jMg)t + (MgtJMy 9 )t Je 


: : . a : 3.) —BG 
+ [(M,-jMy) + (Mo-jM_)t + (Mo-jMg)t +(Mg-JM, Q)t Je ; 


Ce) 
20 


Combining terms Equation (29) becomes: 
-5t -5t 
Kyo + 2M,e cos(Dt) + 2Moe alge Odie) 
a 
a" at 


(2M+2M.t+2M_t +2Mgt>)e ems (De) 


y(t) 


a+ 


aor 


i: (2M +2Mpt+2Mgt~+2M, \t>)e = sin( Dee (30) 


In this case as in the previous two cases considered, gnhe 
steady-state solution is given by Equation (20). In Table 
VI, the reversion solution is compared to a machine solution 
for ao = 1. Ose = 4.0, h = 2.0 and c = WOR in Tableman. 
the solutions are compared for the non-linear coefficient 
equal to ten.. A comparison of the results Iassted in @ables 
VI and VII with the results obtained in the previous cases 
for the same values of we 


reversion method yields a less accurate transient solution 


heand c, indi@aces that the 


when the damping coefficient is small. In this case the 
linear solution is seen to be a better approximation over a 


mortion of the transient solution. 


B. oOLNUSOIDAL FORCING FUNCTION 
Substitution of a constant-amplitude sinusoidal forcing 
function, c Sin(we), into the first two temmeror Equar wan 


(15) and utilization of the trigonometriemidentity [Rema 
ae , oe i 
sain ~ Cw) 7sin (ut ) - 7sin( 3ut ) (31) 


yields: 


eal 


Js ou bnleh 8/1) 


a=1.0, es Oe howe 0 , te=1. 6 
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TABLE VII 


a=1.0, wo=4.0, h=10.0, e=1,0 





Time 


aROY: 


Or: 


SS WN 


BR eS 


Se oa FF HO Be UI UU 


SP Ul oo UI 6S 


U7 


= “oi 


Linear 
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Y(s) = —_So ? 
(s°+u>) (StA)GerE) Oe tu) (shh) GetB) 


pherw 


oe 2) 
<8) hal 9 a = 
As in the case of the constant forcing function three solu- 
tions are considered depending on the nature of the radical 
insfouetioms (is) and (14). 

1. Radical Greater than Zero 


The inverse Laplace transform of Equation (32) is: 


sin(wt) +N et Noe e+ M cos(wr) + 


y(t) =M,,cos(wt) +™M 1 13 


12 


M,ySin(wt) + (N2+N,t+N 6° +N-t- jen 7 


5 


(Ny + Net + Ngt® rn Nj te” + Mj ,cos(3ut) + M,¢sin(3wt) + 


—-Bt 


Q 3, -At 2 3 
(N.. +N,,.t+N..t +N, yt Je + (Ny ot Nit +N i +Nigt Je 


ala Ne as, ay 


(Big) 
The solutions represented by Equation (33) and the first 


four terms of Equation (33) were tested using a = 8.0, 
wi=9.0, c=1.0, h=20.0 and w=1.0. The solutions are compared 
tO a machine solution in Table VIII. In Table IX the solu- 
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to fifty. Table X is a comparison of the solutons for the 
iecaguency Gmmume forcing Temetion increased to five gage 
non-linear coefficient equal to fifty. Inspection of the 
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increased, the accuracy of the solution decreases and as 
the frequency of the forcing function is increased, a 
more accurate approximation is obtained. The effect of the 
magnitude,of the» amplitude. of the forcing function wasmm@ot 
illustrated; however from an Sean nee of the coefficients 
of Equation (33) it can be seen that as c is increased it 
would have a more detrimental effect on the accuracy of 
the solution than corresponding increases of the non-linear 
coefficient. An examination of the solution for the non= 
linear coefficient equal to zero, represented by the first 
four terms of Equation (33), reveals that it is overall as 
accurate a solution as the solution utilizing all the terms 
Of stquatiom (33): 

Oo. Radical Equal to Zero 


Equation (32) becomes: 


( ae phemw pho?w 4) 
(s“¥w")(s#A)* (s?+w2) (s4a)° (s74+9w-) (sta) ® 


The inverse Laplace transform of Equation (34) yields: 


y(t) = M,7cos(wt) +M, gsin(wt) +(Nyg+Nogtie + M, gcos (wt) + 

My gsin Cut )+(N55#N5 ott, ot °4N 5 t54N ot AN pt°4N, tO HN gt eae 
Mp, cos (3ut)4My,sin(3ut)+(N5otN a ttNq,t°#Na5t +N at. + 
Nyt? +Nzot° +Nact!)eA¥ (35) 
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Table XI is a comparison of the solutions for the parameter 
values a=6.0, wi=9.0, h=20.0, c=1.0 and w=1.0. A comparison 
of errors listed in Table XI with the errors listed in 

Table VIII indicates that a decrease in the magnitude of 

the damping coefficient results in a less accurate approxi- 
mate solution. The linear solution represented by the 

first three terms of Equation (35) is again seen to be 
overall as accurate an approximation when compared to the 


solution represented by Equation (35). 


3. Radical Imaginary 


Utilizing the relationships expressed in Equations 
(25) and (26) and combining terms with complex conjugate 
coefficients, the inverse Laplace transform of Equation (32) 
DECOMe s ; 


Ot Oy 


watt) = M5 3cos(wt) +M5,sin(wt) +N - cos (Dt )+Na,¢ é sin (Dt) + 
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ae-~* Ui) 96o2n GW) 39 HO Ng oO eC Cos 
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(Ny gtNy yt+Ny et +Ny ct? )e e sin (Dt det] seme Sail ea, 


ay 
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Mo g8in(3ut )+(Ny7+Ny pt+Ny gt +N nie - cos(Dt) + 


> 
ee 


t+N t°+N t> de - sin@inlae (36) 
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Using the values a=1.0, wo=9.0, h=20.0, =i O sand = ise 

the solutions are compared in Table XII. A comparison of 
Table XII with Tables VIII and XI further illustrates the 
oie tacatl.effiectk..of dkhae dampizae coeffichent on, the accumacy, 


of the solution. Inthe three tables the parameter values are 
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equal except for the magnitmes of the damping coeil i ican 
As a was decreased from eight in Table VIII to six in 
Table XI and finally to one in Table XII, the magnitude 
of the error increased correspondingly. The linear solu- 
tion illustrated in Table XII is seen to be an overall 


better approximation than wsthesreversionysotution. 


ee 


VI. MODIFIED SOLUTION 


The reversion method previously examined was limited in 
application to a region dm Which the solutiong@dad wie: 
diverge. The reversion method can be made to always con- 
verge by modifying the solution to be correct in the steady 
state for any parameter values. The steady state for a 
eenstant forcing function can be @bBTaimed by solving wae 


enore equation: 


wey - hy? = Ce (37) 


Consider a truncated form of Equation (23): 
At -Bt 


y(t) = Koy tK5,e + K55 


-AT 
e + Ky +(K5y +Kyote + 


(Kop tKygt)e (38) 


Since the coefficients K Boek have the common factor 


a5 29 
he, see Appendix, Equation (38) can be rewritten as: 


3 
% -At -Bt he ewe : -AT 
y(t) = Rog + Ko4e + Ko5€ See + He (K5y + Ko5t)e 2 
0 


ne7(Kyg +Kygt)e (39) 


Let SS be equal to the steady state obtained from the solu- 
tion of Equation (37). If the steady-state portion of 


Equation (39) is set equal to SS, then: 


3 
ager 5-4 (40) 


ae 


Equation (40) can be rewritten as: 
eo Sen eS es) (41) 
0 20 


Substitus@tonmeimigquation” = (4") amto Equation C39) yrerae 


a -At -Bt Qi ~At 
y(t) = SS+K,,e + Kj5e + (Wo) (Ky 9-SS)(K5y + Koet)e + 


-Bt (42) 


(we) 


(K5q -SS) (Kg + Kygte 
Equation (42) is the modified solution. Using the values 
a=6.0, we=H.0, h=1.0 and c=1.0, the modified solution is 
compared to a machine solution in Table XIII. The results 
indicate a relatively accurate approximate solution was 
obtained. In Tables XIV and XV the solutions are compared 
for the non-linear coefficient respectively equal to ten 

and fifty. The tables reveal that the error increased as h 
was increased; however, the modified solution is seen to 

be still usable as an approximation for the non-linear 
coefficient equal to fifty. In Table XVI the solutions are 
compared for the forcing function equal to ten. The magni- 
tude of the forcing Pa ee is seen to be a more critical 
factor on the accuracy of the approximation than correspond- 
ing increases of the non-linear coefficient. A comparison 
of the modified solution with the reversion solution for 
Similar parameter values indicates that the modified solution 
is a more accurate approximation for large values of h and c. 
For small values of h and c, the modified and reversions 
solutions are of the same relative accuracy, as revealed by 


a comparison of Table XIII with Table IV. 
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VII. CONCLUSIONS 


A method for obtaining an approximate solution for a 
specific non-linear differential equation has been derived 
utilizing reversion-of-series and non-linear transform 
techniques. The usefulness of the solution was found to 
be dependent upon the nature of the forcing function and 
the parameter values specified. Also, it was found that 
useful approximation could be obtained for particular par- 
ameter values by simply considering as a solution that 
portion of the reversion approximation which was independent 
@f the non-linear coefficient. For the case of the constant 
forcing function the linear approximation was found to be 
less aecurate than the reversion solution. In the case 
of the sinusoidal forcing function the linear approximation 
was found to be as accurate as the reversion solution when 
the damping was significant and a more accurate approxima- 
tion when the damping coefficient was relatively small. 

For a constant forcing function it was found that the 
steady state could be specified and utilized to modify the 
Peversion solution. The modified solution was foumd to be 
useful as an approximation over a relatively wide range of 
parameter values and gave results that were superior to The 
linear and reversion approximations. 

The reversion method was used to solve a particular non- 


linear differential equation; however, the development of the 
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solution indicates the possibility of applying this method 
to solve non-linear differential equations Similar in forum 


to the equation considered. 
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